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Abstract 



We study the maximally supersymmetric Kondo model obtained by adding a fermionic 
impurity to M = 4 supersymmetric Yang-Mills theory. While the original Kondo prob- 
lem describes a defect interacting with a free Fermi liquid of itinerant electrons, here the 
ambient theory is an interacting CFT, and this introduces qualitatively new features into 
the system. The model arises in string theory by considering the intersection of a stack 
of M D5-branes with a stack of N D3-branes, at a point in the D3 worldvolume. We 
analyze the theory holographically, and propose a dictionary between the Kondo problem 
and antisymmetric Wilson loops in Af = 4 SYM. We perform an explicit calculation of 
the D5 fluctuations in the D3 geometry and determine the spectrum of defect operators. 
This establishes the stability of the Kondo fixed point together with its basic thermody- 
namic properties. Known supergravity solutions for Wilson loops allow us to go beyond 
the probe approximation: the D5s disappear and are replaced by three-form flux piercing 
a new topologically non-trivial S 3 in the corrected geometry. This describes the Kondo 
model in terms of a geometric transition. A dual matrix model reflects the basic properties 
of the corrected gravity solution in its eigenvalue distribution. 
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1 Introduction 

The Kondo model has been a rich subject of study for almost 50 years. In a certain sense, 
it was the first system to exhibit flow from an asymptotically free UV fixed point to a 
strongly coupled IR regime. It has also been a playground for early renormalization group 
studies, exact conformal field theory techniques, and large N methods in field theory. Its 
generalization to the Kondo lattice system yields models which may be of great relevance 
in understanding non-Fermi liquid behavior in the heavy fermion metals. For reviews with 
extensive references, see [TJ [2] . 

A supersymmetric variant of the Kondo system, which may be useful in attempts to apply 
AdS/CFT to condensed matter systems, was described in [3j (see also pfl EJ [6] for related 
work). The analogue of the Kondo defect spin arises at an intersection of M D5-branes 
with a stack of N D3-branes. At large 't Hooft coupling for the D3-brane field theory, the 
D5s appear as probes in an AdSc, x S 5 geometry, wrapping an AdS2 x S 4 submanifold. They 
describe fundamental defects interacting with the ambient SU(N) M = 4 supersymmetric 
gauge theory (which provides the analogue of the "itinerant electrons" in the Kondo model). 
Unlike the original Kondo model, in our case the ambient field theory is already at a 
nontrivial interacting fixed point, giving rise to qualitatively new features that we will 
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explore here. 



In this paper we continue the study of this maximally supersymmetric Kondo system, 
extending the results of j3] in several directions. We will make explicit the correspondence 
between Kondo impurities and half BPS Wilson loops of = 4 SYM, and this will 
allow us to extract new results at strong coupling. Some of the main questions about the 
Kondo model involve the physics reflected in the defect entropy specific heat, and magnetic 
susceptibility. In the latter two cases, the first non-trivial contribution involves knowledge 
of the leading irrelevant operator (invariant under any relevant global symmetries, like spin 
and channel symmetries) at the IR fixed point [TJ |2]. So to determine the basic physics, 
one must find the spectrum of operators in the defect theory, which maps in the gravity 
regime to a classification of KK modes on the probe D5s. The complete spectrum and 
map to the dual field theory has not yet been determined; we complete this task here. 
Our results reveal that the superconformal algebra which governs the fixed point places 
powerful constraints on the defect thermodynamics (and can even strongly limit the effects 
of inter-defect interactions, in models with multiple D5 branes). 

While gauge/gravity duality immediately sums the effects of the large N gauge theory on 
the defects (through the details of the geometry wrapped by the D5-brane probes), the 
leading approximation does not include the backreaction of the defect flavors on the M = 4 
supersymmetric field theory. To determine this at strong 't Hooft coupling, one must find a 
corrected solution to the supergravity equations, including the non-trivial effects of the D5- 
brane stress-energy sources instead of simply treating the D5s as probes. These solutions 
are already known, from earlier work of D'Hoker, Estes and Gutperle [7]. We describe them 
explicitly here and apply their results to the analysis of our supersymmetric Kondo model, 
now including backreaction. The D5-branes induce a geometric transition, changing the 
global topology of the solution, and replacing themselves with three-form flux piercing a 
new non-trivial 5* 3 in the corrected geometry. We also discuss how the physics of these 
backreacted solutions is captured by a matrix model [8], which may be useful in carrying 
out further field theory computations. 

In addition to their intrinsic interest, the questions we investigate here are well motivated 
because understanding the physics of the single-defect model is the logical starting point for 
a more ambitious program of studying gravity duals to lattices of such defects, analogous 
to Kondo lattice systems. Kondo lattice models of defect spins interacting with a bulk 
Fermi liquid are thought to give rise to the fascinating non-Fermi liquid phases in the 
heavy fermion systems (which in turn may be closely related to the "strange metal" phase 
of the cuprate high temperature superconductors). The papers [9] initiated a holographic 
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approach to non-Fermi liquids using special properties of the AdS2 geometry that emerges 
in the near-horizon limit of the Reissner-Nordstrom black brane. Lattice models of defect 
flavors provide an alternative method of obtaining the geometries that can give rise to 
non- Fermi liquid behavior jH ITU], with some important differences from the Reissner- 
Nordstrom case; the AdS2 geometries now arise on the probe world-volumes themselves. 
General discussions of why Kondo lattice models could be relevant to the problem of strange 
metals from a field-theoretic perspective appear in [UJ [12] . 

The organization of this paper is as follows. In §2, we describe basic field theoretic aspects 
of the Kondo model. In §3, we review and extend the work of [3J, discussing the simplest 
results which are visible from a probe approximation. In §4, we turn to a determination 
of the full spectrum of operators in the defect theory; the new defect-localized operators 
are found by performing a KK reduction on the wrapped probe branes. The extended 
superconformal algebra present in this theory plays a powerful role in organizing the results. 
With this list of operators in hand, in §5 we turn to its implications for the basic questions 
in Kondo physics - questions of the IR stability of the Kondo fixed point, and the leading 
non-trivial contributions to defect specific heat and magnetic susceptibility. Finally, in 
§6, we describe the supergravity solutions which incorporate backreaction of the defects 
on the bulk, as well as a dual matrix model which reproduces the main features of the 
supergravity solutions in its eigenvalue distribution. §7 contains our concluding remarks 
about promising future directions. Some unpalatable details of the computation of the 
spectrum in §4 have been relegated to the appendix. 



2 Field theory analysis 



This first section is devoted to a field theory analysis of the Kondo problem and its maxi- 
mally supersymmetric version. After reviewing the original multichannel Kondo model, we 
describe generalizations based on interacting conformal field theories (CFTs), with a focus 
on the supersymmetric case that arises from M = 4 SYM coupled to a defect fermion. The 
quantum theory is studied in two complimentary ways. First we integrate out the interac- 
tions with the ambient theory and obtain an effective quantum mechanics description of the 
impurity at low energies. This provides a simple way of characterizing the RG evolution of 
the Kondo coupling and of evaluating correlation functions of the defect fermions. On the 
other hand, we may also integrate out the defect fermions and derive a new ambient theory 
where the impurity has been replaced by certain boundary conditions on the bulk fields. 
Their equivalence will reveal a very useful correspondence between Kondo impurities and 
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Wilson loops. 



2.1 Multichannel Kondo model 



The multichannel Kondo model is described by the Hamiltonian 

P> a pp' ia/3 

where ipp ia is the creation operator for a conduction electron with channel index i = 



1, . . . , K, SU(N) spin a, and energy 

f ® = |l- fF . (2.2, 

The interaction is antiferromagnetic, J > 0. These "ambient" electrons interact with a 
spin S a/ 3 localized at the origin and transforming under a representation R of SU(N). The 
Fourier transform of this delta-function interaction gives rise to the sum over momenta in 



the second term of (2.1 ). The SU (N) generalizes spin, and this extension will be important 



for us because the system that we will study has N = K (large). 

Our analysis will focus on the antisymmetric representation with k indices, in which case 
the impurity spin S can be described in terms of auxiliary fermions Xa with total fermion 
number k: 

k A 

Sa^ = xlx? ~ T7<W > Z2 X ^ Xa = k ■ ( 2 ' 3 ) 



a=l 



The original version of the model corresponds to one channel (i.e. K — 1) with SU(2) 
spin (N = 2), interacting with a single impurity fermion. Integrating out the conduction 
electrons up to some cutoff of order of the Fermi momentum shows that at low temperatures 



the Kondo coupling flows to J — > oo. This means that the Hamiltonian (2.1) is minimized 
by binding one electron with the impurity to form a spin singlet (recall that J > 0). 
Therefore, the impurity spin is completely screened and disappears from the low energy 
theory. We are left with a theory of free electrons, but we have to impose the boundary 
condition ip = at the origin. This phenomenon will have an analog in the gravity dual, in 
terms of a geometric transition where the D5 branes that give rise to the impurity disappear 
and are replaced by their flux. 

In the multichannel K > 2 case, the ground state of the strong coupling point J — > oo 
is a Young tableau with K — 1 columns and N — k lines; this has larger dimension than 
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the degeneracy at zero coupling, given by N\/k\(N — k)\ (see e.g. [13] ). Therefore, a 
RG flow from weak to strong coupling is not possible, and instead an intermediate fixed 
point is expected. This is known as "overscreening" , and is not described by Fermi-liquid 
theory. Our supersymmetric model will correspond to this regime, and we will study the 
intermediate fixed point using AdS/CFT. 

Let us now summarize some of the results on the Kondo model, obtained using CFT and 
large N techniques. In the approach of Affleck and Ludwig (see j2j [1] for review and 
references), the system is first mapped to a 2d theory by performing an s-wave reduction. 
The ambient CFT is decomposed into spin, flavor and charge degrees of freedom, and is 
described by a Kac-Moody algebra SUk(N) s <g) SU n{K)j <g) U(l) c . At the intermediate 
fixed point, the local impurity 5* disappears, becoming part of the spin current. Thus, its 
correlator is fixed by conformal invariance, 

WWI-i^-^' {2A) 

The whole effect of the impurity is encoded into constraints that determine the number of 
times a given state occurs in the spectrum. According to the "fusion rule" , the spectrum 
is determined by acting on a given primary operator with the spin current in the repre- 
sentation R of the impurity spin. The flavor and charge sectors are not modified. These 
constraints give rise to a conformal field theory with boundary conditions that respect 
conformal invariance -a boundary CFT. 

The approach of integrating out the ambient electrons to obtain an effective theory on the 
defect has also been followed in the literature; see e.g. [13]. The saddle point equations for 
the defect fields can be solved at large N, leading to results consistent with the boundary 
CFT predictions. 

We will then compare these to the predictions from the gravity side. Particularly important 
is the impurity entropy (or g-function), which is defined as 

logg = Simp = liniT^o limy^oo [S(T) - S amhie nt(T)} (2.5) 

where iS am bient is the contribution proportional to the volume V. The g-function is the 
analogue of the c-function for theories with boundary. In the multichannel model, 

k • 7r(Af+l-n) 

sm — — 

<S imp = log J] ^ N Jl ■ ( 2 -6) 

„=1 b111 N+K 
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2.2 Supersymmetric Kondo model 



The Kondo fixed point is naturally analyzed in terms of a 1 + 1-dimensional CFT, obtained 
by performing an s-wave reduction 

ip(p) = rfoip) + higher harmonics (2-7) 

and keeping only the spherically symmetric wavefunctions. Linearizing around the Fermi 
surface and using the doubling trick obtains 

/+oo 
dx ifj[ (id x + J5(x) S) i) L . (2.8) 
-oo 

The model reduces to a 2d CFT coupled to degrees of freedom living on the boundary. A 
natural extension is to couple an interacting CFT to a localized impurity, 

S = S CFT + J 'dtx(idt + O{x = 0)) X , (2.9) 

where O is some operator of the CFT, and the impurity is at the origin x = 0. We refer the 
reader to the review [12] and references [H] for works on this approach in the condensed 
matter literature; in particular the latter papers study Kondo defects coupled to CFTs, 
albeit simpler CFTs than the M = 4 super Yang-Mills theory. 

Our AdS/CFT techniques will be closest to this approach, albeit with additional super- 
symmetry and in the large N limit. We couple M = 4 supersymmetric SU (N) gauge theory 
(the ambient fields) to M (0 + l)-dimensional fermions Xi with total fermionic number k 
(the "defect"). The action reads 

S = S N=i + J dt i [xld tX ) + xl(T A ))(A (t,0) A + n a( p a (t,0) A )x i I \ , (2.10) 

where i, j — 1, . . . , N, I — 1, . . . , M. Here T A are the generators in the adjoint of SU(N), 
and n a is a unit vector in IR 6 . We also need to impose the constraint Yli XiXi = which 
fixes the fermion number. It will be convenient to impose this using a Lagrange multiplier 
M x M matrix Aq, and adding 



S k = Jdt {xKAoYjX'j ~ fc(io)j) , (2.H) 



to the action (2.10). In the brane construction below, A arises from the D5 gauge field. 



In the conventions of (2.10), the kinetic term for the gauge field and scalars is proportional 
to i/gYMi so t ne coupling between the impurity and bulk degrees of freedom is proportional 
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to the gauge coupling Qym- The N = A theory preserves 16 supercharges e; half of these 
are broken by the defect, 

7o n a 7 ae = e, (2.12) 

where 7m are the lOd Dirac matrices. The defect fermions are singlets under the preserved 
supersymmetries. The large amount of supersymmetry, which fixes the Kondo coupling in 
terms of gYM, has important consequences on the quantum theory. We will argue shortly 
that the defect action does not receive quantum corrections; furthermore, having a fixed 
Kondo coupling restricts the spectrum of defect operators, as will be seen in 



Comparing with the multichannel model (2.1), we see that our theory corresponds to 



K = N. This particular case turns out to be quite interesting on its own, because matrix 
model techniques can be applied to the Kondo problem, as we discuss in more detail in £j6j 
On the other hand, the defect fermions x\ have an additional U(M) flavor number that 



is absent from (2.1). This allows for more general (antisymmetric) representations for the 
spin impurity. 

In £j3] we will review how this field theory arises from N D3 branes intersecting M D5 
branes (with worldvolume flux k) at a point. Taking the limit of large N and strong 
coupling will then allow us to solve for many properties of the theory using the gravity 
dual. Before moving in this direction, in the rest of the section we discuss basic properties 
of the quantum theory at weak coupling. 



2.3 Quantum effects from the ambient CFT 



Let us integrate out the ambient fields to obtain an effective + 1 dimensional theory. We 



start from (2.10) at weak coupling and zero temperature, and integrate out A + n a (p a . The 



effective action on the defect will be of the form 

Se ff = J dtt X ld tX t I + 3 ^ I dtdt'D{t-t')x{(t) X ] Am\t')itj{t') + S k , (2.13) 
where the kernel D(t — f) is generated by the scalar and gluon propagators, 



(A (t) A A (t') B + n a n b (j) a (t) A( j) b (t') B ) 



9ym 
4vr 2 



8 AB D(t-t'). 



(2.14) 



A, B are indices in the adjoint, while a, b run over the six scalars of = 4 SYM. Also Sk 



was defined in (2.11) and determines the fermion number. 



To continue, we need the explicit form of these propagators, which at tree level (and in 
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Feynman gauge) is 



(A A Ax)A B (x')) 



9yM $ABS ab 



47T 2 



(X — X 



A2 



9ym Vh^ 



AB 



47T 2 



[X 



X 



A 2 



(2.15) 



where ^ = diag(— 1, 1, 1, 1). Replacing (2.15) into (2.14) reveals something quite striking: 
the contributions from Aq and exactly cancel, giving D(t — if) = 0! In fact, this is also 
true to all orders in the gauge coupling JTOJ El HE] • We conclude that the defect theory 
is not renormalized at zero temperature. This rather special behavior is due to the large 
amount of supersymmetry, which protects the defect action against quantum corrections 
and fixes the Kondo interaction in terms of gyM- 

Next we put the theory at finite temperature T = 1//3 and analytically continue to imag- 
inary time r = it. This breaks supersymmetry and we expect nontrivial quantum correc- 
tions]^] The propagators at finite temperature can be obtained from the zero temperature 
limit ( |2.15 ) by an inversion that maps the line (the time direction) to a circle. The im- 
portant point here is that while the scalar propagator is invariant, the gluon contribution 
is not. The gluon propagator is just changed by a total derivative [T7] or, equivalently, a 
constant contribution to D{r — r') after replacing into the effective action. On dimensional 
grounds, D ~ l//3 2 ]^] In contrast, the ambient fermions in the Kondo model of {2.1 have a 
propagator that scales like 1/r. 



In summary, the defect theory at finite temperature becomes 

S eff = 



xfrx} • / drdT'xl{T) X ^r)x\r') jX ],{r') + S k . 



32/3 2 



(2.16) 



Equipped with this action, we are now ready to calculate the fermion correlator and im- 



purity entropy, following [121 Ej. From (2.11), A plays the role of the chemical potential 
for the defect fermions; let's denote its expectation value by 

1 rf> 



P 



driAo 



(2.17) 



The quartic interaction gives a contribution to the self-energy S oc D(t)G(t), so at large 
N the Schwinger-Dyson equation for the two-point function becomes 

1 



G(u) 



LO — fM 



X 

16/3 



(2.18) 



1 A recent related discussion has appeared in [5], where the behavior of D(t) was conjectured based on 
comparisons with the D5 probe limit. Here we derive D(t) directly from the field theory, thus verifying 
the assumption of [S] and obtaining nontrivial agreement with the gravity prediction. 

2 This may also be calculated from the circular Wilson loop, as explained e.g. in [17] . 
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where the 't Hooft coupling A = Ng^M comes from the quartic interaction and the N- 
component fermion in the loop. The solution is 

l/2> 



G{u) 



A 



16/3 2 



CO — /! — 



A 



16/3 2 



(2.19) 



The result (2.19) determines the basic properties of the system. The chemical potential is 
fixed in terms of the fermion number k by G{r — > 0) = k/N. At large 't Hooft coupling 
this yields [5] 



k 
N 



1 
2 



1 



arcsm /i 



2/3 _ 



(2.20) 



7T 7T 

In particular, this implies that the chemical potential increases with A as jj, oc a/A//3. Now, 
recall that /2 is related to the defect free energy by ft = —dF^eiect/dk; deriving once more 
with respect to the temperature obtains 



imp 



dk 



dfx 
df 



Using the result (2.20) we thus arrive to the impurity entropy 

Ny/X 



-"imp 



3tt 



1-fi 



2\3/2 



(2.21) 



(2.22) 



The same result will be obtained from the gravity side, with \x = cos 6k being related to the 
angle that defines the embedding of the D5 branes inside the S 5 . Notice that the impurity 
entropy increases with the 't Hooft coupling -an effect that is absent in the original Kondo 
model. From the field theory calculation, this is due to the fact that the chemical potential 
increases like VX 



2.4 Impurities as Wilson loops 

We can also derive an effective theory for the bulk fields by integrating out the defect 
fermions, which can be done exactly. This produces a Wilson loop insertion, 

W R = Tt r P exp (i J dt(A + ny a )^ (2.23) 

in a representation R of the gauge group. This establishes a correspondence between the 
supersymmetric Kondo model and Wilson loops in = 4 SYM. In the rest of the paper we 
will develop this dictionary in detail, and use exact results on Wilson loops to understand 
the strong coupling behavior of the Kondo modelj^] 

3 A similar relation between Kondo-CFT problems and Wilson loops was discussed in the condensed 
matter literature in e.g. [T9] , 
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Let us review how the Wilson loop is generated in the simplest case M — 1; we follow |15j . 
It is convenient to choose a gauge where A + n a (p a has constant eigenvalues mi, . . . , m^. 
The equation of motion for the defect fermions is 

(td t + mi ) X i = , i = l,...,N. (2.24) 

We have a Fock space of N fermions, each with energy mj, and in the partition function we 
only need to sum over states with k fermions. The partition function of the defect fermions 
then becomes 

Aiefect = 2^ e • • ' e ' (2.25) 

ii<-»a<— <*fc 



where an infrared regulator is assumed. Reintroducing the combination A + n a (j) a , (2.25) 
is recognized as the trace of the Wilson line in the k-th antisymmetric representation of 
SU(N), 

e^ dtm ^ . . . e ^ dtm ^ = Tr Ak P exp (i j dt(A + n a (f) a )) . (2.26) 

This shows that integrating out the degrees of freedom that live on the defect produces a 
supersymmetric Wilson loop insertion in the M = 4 theory. 

In the case M > 1 we have NM fermions with total fermion number k. This allows for 
general representations for the Wilson line, where the Young tableaux has at most iV rows 
and M columns, and the number of boxes is determined by k. (For M = 1 this simply 
gives k rows). Now the matrix structure of A needs to be specified as well; choosing 

A = dmg{Q 1 ,...,Q M ), (2.27) 

obtains an equation of motion 

^O t + mi + n I ) x I i =0. (2.28) 

So there are NM fermions, with energy m, + Qi and total fermion number k. The chem- 
ical potential Qj determines the number of boxes in the J-th column as the value of the 
propagator (x/( r )x J (0)) as r — > 0. 



3 D3-D5 system in the probe approximation 

Having understood the basic properties of the field theory with a quantum impurity at weak 
coupling, in what follows we will study the Kondo model at large 't Hooft coupling using 
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AdS/CFT. We begin by reviewing the realization of quantum impurities using intersecting 
D3 and D5 branes. Next we study the model at large 't Hooft coupling and large N, 
but keeping the number of D5 branes M fixed. As explained in [3], this is equivalent to 
having D5 branes embedded in AclSc, x S 5 with worldvolume flux. From the field theory 
perspective, this is a quenched approximation where loops of defect fields are neglected. 
After reviewing the gravity description, we calculate the g-function and compare with the 
nonsupersymmetric Kondo model of §2.1 



3.1 D-brane description of the Kondo model 

The brane configuration is given bjj^] 








1 


2 


3 


4 


5 


6 


7 


8 9 


N 


D3 


X 


X 


X 


X 












M 


Dh 


X 








X 


X 


X 


X 


X 


k 


Fl 


X 
















X 



This configuration preserves 8 supercharges. The k fundamental strings between the D3 
and D5 branes lead to k fermions localized on the D3/D5 intersection. 



The full open string action is given by 

S = Sd3 + Sdc, + Sdefect 



(3-1) 



where Sd3 is the action for the N D3-branes (namely Af = 4 SYM), Sd5 describes the 6d 
worldvolume theory on the M D5-branes, and Sdefect is the action for the 3-5 strings at the 
intersection, 



S, 



defect 



dt 



rfAx\ + xl (A (t, 0)5 + n a <p a (t, 0)}) X \ + xl(Ao)iti ~ M4>tfl (3-2) 



where i, j = 1, . . . , N, I = 1, . . . , M. The fermion x is a bifundamental of SU(N) x SU(M), 
so it couples to the D3 and D5 gauge fields A and A . Supersymmetry requires that it 
couple also to one linear combination n a (f) a of the scalars of the M = 4 theory; n a is a unit 
vector in IR 6 . Aq also acts as a Lagrange multiplier fixing the fermion number xlx l j = k5j. 

In the decoupling limit for the D3 gauge theory, the 5-5 strings become nondynamical. Then 
we are left with 4d Af = 4 SYM with gauge group SU(N), coupled to a + 1-dimensional 



We refer the reader to [H [15] for detailed descriptions of this system. 
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fermionic system (3.2) with flavor symmetry U(M). This reproduces the supersymmetric 
Kondo model of 52.2. 



Let us now discuss the symmetries preserved by the defect. A line (x°, 0, 0, 0) preserves 
rotations along the line, time translations, dilatations and special conformal transfor- 
mations along the time direction, 6x° = b°x 2 . These generators (Jjj, P , D, K ) give 
SO(3) x SL(2,WL); this is the subgroup of the original 50(4, 2) conformal group that 
is left unbroken by a line. On the other hand, the £0(6) R-symmetry that acts on the 
Af = 4 SYM scalar fields is broken to SO (5) by the coupling x n a(p a X- Thus, the bosonic 
symmetries preserved by the Wilson loop are SL(2, R) x 5*0(3) x SO(5). The defect also 
breaks half of the supersymmetries of the ambient gauge theory. These combine with the 
bosonic symmetries to give the supergroup OSp(A*\A) [15J. This group plays an important 
role in understanding the dynamics at strong coupling. The operators of lowest dimension 
correspond generically to short multiplets of OSp(4*\A), whose dimensions are protected 
and determined by group theory. We calculate the spectrum of excitations explicitly in £|4| 
and use these results in ^5] to study the stability and thermodynamics of the defect theory. 

The dynamics of the system is controlled by the 't Hooft coupling A = 9ym^ anc ^ M/N. 
At weak coupling A < 1 we can perform perturbative field theory calculations, and the 
results of £|2] apply. On the other hand, at large 't Hooft coupling there is dual gravity 
description in terms of type IIB supergravity with 3- and 5-form fluxes. When M/N <C 1 
the D5 branes can be treated as probes in the AdS§ x S 5 background generated by the 
D3 branes. For 1 <C g YM M <C 9ym^ the backreaction of the 5-branes should be taken 
into account as well. We now focus on the probe limit, postponing the analysis of the fully 
backreacted solution to £j6j where we will make contact with the results of [7]. 

3.2 Gravity side 

Let us review the gravity description of Wilson loops in the probe limit, where M/N correc- 
tions are ignored, and apply it to the Kondo model. As found in [3J, in this approximation 
the quantum impurity is described holographically by D5 branes embedded in AdS§ x S 5 , 
with worldvolume flux. In field-theoretic terms, this gives a solution to the Kondo model 
in the quenched approximation, where loop effects of the defect fermions are ignored. Fur- 
thermore, the ambient conformal field theory is very strongly coupled, and its quantum 
effects are resummed into a gravitational dual. This approximation will allow us to find 
important quantities such as the defect entropy and thermodynamics controlled by the 
leading irrelevant operators. 
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In order to realize the defect symmetries transparently, it is useful to slice AdS$ x S 5 by 

ds 2 = R 2 (du 2 + cosh 2 u ds 2 AdS2 + sinh 2 u dVt 2 2 + d6 2 + sin 2 9 dfi 2 ) , (3.3) 

with u > and < 9 < tt. The potential Ca that gives rise to the five form flux can be 
written as [8] 



C 4 = R 4 



■- + - sinh 4 M ) eVeV + (-9-sm29+- sin 49 ) e 6 e 7 e 8 e 9 
2 8/ \2 8 



(3.4) 



where (e°, e 1 ) are the vielbeins for ds\ dS2 , (e 4 , e 5 ) are the vielbeins for dfl 2 ,, and (e 6 , e 7 , e 8 , e 9 ) 
are those for dQ 2 . 

The D5-brane worldvolume is AdS2 x S* 4 , given by the embedding conditions 



M = , 9 = 9 k . (3.5) 
As we review shortly, this angle is related to the number of defect fermions k by 

/»■ = — I H, - - sin 29,.. ) . (3.6) 



N ( 1 

— fc - - sin 26> fe 

7T V Z 



This configuration was studied originally by 



The bosonic symmetries SL(2, 



x 



5*0(3) x SO (5) are explicitly realized in the way of writing the metric (3.3). 
Ignoring nonabelian interactions, the D5 branes are described by the DBI action 

S 5 = -T 5 J d^i v 7 - det(G + F) + T 5 J F A C 4 , (3.7) 



with worldvolume metric 



ds 2 m = R 2 (ds 2 AdS2 + sin 2 9 k dVt\) . 



(3.J 



The background value (|3.4|) induces a tadpole for F due to the Chern-Simons interaction. 

(3.9) 



As a result, there is nonzero worldvolume flux given by 



F = cos 9 k e° A e 1 , 
and the fundamental string charge k is given by SSd 5 /5F i. 



3.3 Defect free energy and entropy 

As a first step towards solving the supersymmetric Kondo model, we now evaluate the 
boundary free energy and entropy (g-function) in the probe approximation, and compare 
with the multichannel Kondo model of 32.11 
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The defect free energy, as a function of temperature T, can be computed using the DBI 
action immersed in an AdS black hole. This calculation is easier to perform in the original 
Poincare slicing of AdS, where the black hole metric obtains the familiar form 

„ , * . n dr 2 ^ 2 i - 3 



-f(r)dt 2 + — + — ( dx i ) + R2 ( dd2 + *™ 2 6 dQ 



4) , 



(3.10) 



,i=l 



(3.11) 



where 

M = W (l " ^ 

The temperature of the field theory is related to the location of the horizon by T = r + / ttR 2 . 
The 5-form flux is 



F 5 = -(1 + *)-^4 (dt A dr A dx\ A dx 2 A dx 3 ) . 
In these coordinates, the D5 branes are extended along r, t and the S A . 



(3.12) 



After regularizing by subtracting the Euclidean action of the analogous D5 embedded into 
the pure AdS§ x S 5 spacetime, the D5 free energy at leading order in 1/N is given by 



defect 



sin 3 9 h 
3tt 



NT 



(3.13) 



where A is the 't Hooft coupling A = gy M N. Ignoring nonabelian interactions, for M D5 
branes the total free energy is simply M times the above result. Recalling the field theory 



definition (2.5) of impurity entropy and that S = —dF/dT, from (3.13) we obtain 



log g = S imp = yf\ 



sin 3 9 k 
3tt 



-MN. 



(3.14) 



in agreement with the field theory result (2.22). We note that this coincides with the 



expectation value of the circular Wilson loop, up to a factor of 2 from mapping the Polyakov 
to the circular loop [5]. 

This expression, which is valid at large 't Hooft coupling and when N/M <C 1, is indepen- 
dent of temperature. The result applies to any nonzero (even very small) temperatures and 
is dictated by conformal invariance. On the other hand, CFT techniques also determine 



the defect entropy of the overscreened Kondo model, given in (2.6). At large N, this can 
be approximated by [13] 



-"imp 



-MN 



71 



f 



71 



1 + 



K/NJ 



f 



71 



:(l-k/N))-f 



7T 



:k/N 



(3.15) 



l + K/N" 1 7 J \l + K/N 

where / (x) = f Q du log sin u and recall that K is the number of channels of the ambient 
electrons. 



16 



A crucial difference between both expressions is that (3.14) increases with the 't Hooft 



coupling, which fixes the strength of interactions in the ambient CFT. However, the multi- 



channel Kondo model of [2^ is based on a gas of noninteracting electrons, and this effect 
is absent. Physically, the increase of logg with a/A in our supersymmetric model is due 
to the fact that the chemical potential for the defect fermions is proportional to This 
may be understood directly in field theory terms, as we found in §2.3[ It implies that the 
g-function increases under a marginal deformation of the ambient theory by the operator 
F^ v that changes qym ~a property that could also be understood using the approach of 



To continue comparing the two answers, note that both g-functions become largest when 
k/N = 1/2, corresponding to having particle-hole symmetry x ^ X- Normalizing the 
g-functions so that their value is one at the maximum, Figure [T] shows log g as a function 
of k/N by numerically solving (3.14) (dotted line) and from (3.15) for the values K/N = 1 
(blue curve) and K/N = 1/10 (red curve). In particular, this reveals that the two entropies 
are most similar when K = N, namely when the ambient electrons are N x N matrices 
-as we would expect when comparing to the Af = 4 model. 



Log 



gmaxy 




1.0 N 



Figure 1: Impurity entropy as a function of k/N for the supersymmetric model (dotted curve) and 
nonsupersymmetric multichannel model with number of channels K/N = 1 (blue) and K/N = 0.1 
(red). 



It is also useful to compare the two behaviors in the regime k/N <C 1. Expanding (3.14) 
in powers of k/N obtains 



log# = -kMV\ 



1-1 

10 



.~2N 



2/3 



3 
280 



3irk\ 



4/3 



2N 



+ 



(3.16) 
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This intriguing nonanalytic dependence was already noticed by [22} [8] . Interestingly, the 
nonsupersymmetric CFT result also shows a nonanalytic dependence, albeit a logarithmic 
one: 

' TX k \ 7T 2 / k 

2NJ ~ 36 ViV 



log g = fcM 



1 - log 



+ • 



(3.17) 



This difference is related to the fact that one ambient CFT is strongly coupled SYM 
while the other one is a system of free fermions. It would be interesting to have a better 
understanding of these behaviors. 

Lastly, we mention that there is a discontinuity at precisely zero temperature, since any 
nonzero value is equivalent to a large temperature phase. The zero-temperature entropy 
can be computed using the fact that the configuration is supersymmetric and hence one 
just needs to count the number of supersymmetric ground states with the appropriate 
defect fermion number. As in j3], the answer is: 



S imp [T = 0] = log 



k\(N-k)\ 



(3.18) 



This result is exact, and is related to the fact that the expectation value of the straight 
Wilson line is protected by supersymmetry from quantum corrections. 



4 Spectrum of defect operators and holographic excitations 

The M = 4 theory with a defect has gauge-invariant operators localized on the defect, 
which contain the defect fermion \. The dimensions of these operators determine impor- 
tant properties of the Kondo model, such as the stability of the fixed point and defect 
thermodynamics. This section is devoted to the analysis of defect operators and their 
spectrum. After reviewing their basic properties using the OSp(4*\4) algebra, we present 
an explicit calculation of the anomalous dimensions in terms of KK fluctuations of the D5 
branes in the gravity dual. These results will be used in §5ko study the stability and ther- 
modynamic aspects of the Kondo fixed point. The basic properties of the KK spectrum 
have recently been conjectured in [B] using group theory analysis. Here we build upon their 
results, performing an explicit calculation of the D5 fluctuations. This will also show what 
representations of the algebra are physically realized. 
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4.1 Defect operators 



Starting from the action ( 3.1[ ) and (3.2), at large iV the defect operators of lowest dimen- 



sions are single traces of SU(N) indices and contain one pair of defect fermions, 

(o^)/ = j2xK<P ai ...r n )Ui (4-i) 

where <p a are the TV = 4 scalars. These operators transform in the adjoint of the flavor 
U(M) group. Focusing on the spectrum at large N, the operators of lowest dimension 
arise from small multiplets of the 05p(4*|4) superconformal algebra preserved by the de- 
fect. These correspond to chiral primaries, whose dimensions are protected from quantum 
corrections. The even subgroup is SX(2,R) x 50(3) x 5*0(5), and states are classified by 
the corresponding quantum numbers (h, j; mi, 1712), where h is the SL(2, R) dimension, j 
is the 50(3) spin, and (mi,m 2 ) are the 50(5) Dynkin labels of the state. Lowest weight 
representations for this supergroup were analyzed in 



The operator of lowest dimension arises from an ultrashort supermultiplet of 05p(4* |4). 
The full supermultiplet is 

(1, 0; 0, 1) © (3/2, 1/2; 1, 0) © (2, 1; 0, 0) . (4.2) 

The first state has dimension 1, is an 5*0(3) singlet and transforms as a vector of 50(5). 
The next state is a spacetime fermion, and transforms in the spinor representations of 
50(3) and 50(5). The last state, of dimension 2, is an 50(3) triplet and 50(5) singlet. 
For our purposes, we will also need short multiplets with lowest weight state of the form 
(/,0) under SL(2,R) x 50(3); the multiplet content for / > 2 is 

(/, 0; 0, /)©(/ + 1/2, 1/2; 1, / - 1) © (/ + 1, 1; 0, / - 1) © 

(/ + 1, 0; 2, / - 2) © (/ + 3/2, 1/2; 1, / - 2) © (/ + 2, 0; 0, / - 2) . (4.3) 



4.2 D5 fluctuations and KK reduction 

Having understood the multiplet structure of the superconformal algebra, the next step is 
to find which of these operators are physically realized. These defect operators are dual to 
fluctuations of the D5 worldvolume fields in the AdS 5 x 5 5 background. So our strategy 
will be to calculate the spectrum of D5 fields explicitly, and use these results to determine 
the anomalous dimensions of defect operators via 

h scalar = - + -Vd 2 + 4:m 2 
2 2 
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h vector = ^ + ^(d-2) 2 + 4m 2 

h spinor = ^ + \m\, (4.4) 

with d — 1, corresponding to an effective theory on AdS2- Our analysis will focus on bosonic 
fluctuations, since their fermionic parterns may be obtained using the superconformal 
generators. 



The spectrum of fluctuations is obtained by linearizing the D5 action (3.7) around the 



embedding (3.5). The worldvolume metric is (3.8), and there is nonzero worldvolume flux 
(3.9). The relevant fields for us are the 6d gauge fields A^, and the fluctuations 5u and 



59 around (3.5). The AdS 2 gauge field and 59 were studied by [2D]. The full analysis is 



presented in the Appendix, and here we summarize the main results. 
We denote the AdS 2 and S 4 metrics by 

ds\ dS2 = ~ dt2 + ^ , dnl = gpdCde . (4.5) 

The 6d worldvolume coordinates are chosen to be £° = t, = r and (the angles on S A ). 
The relevant fluctuations are 

dt A dv 

5{ds 2 ) = 5u 2 ds 2 AdS2 + (sm29 k 59+cos29 k 59 2 )dnl+(d5u) 2 + ( y d59) 2 , F = cos9 k — + /, 

(4.6) 

where / is a 2-form along 6d space. At this stage, the fields 5u, 59 and / are functions of 
the 6d coordinates £ M . 



The DBI action up to quadratic order is found to be (see (A. 8)) 

,3 



S (2) = T 5 J d & i^) \^~^\ [r 2 ^) 2 - (dM 2 - fyMa - 25u 2 + A59 2 ] 

+ ^ (r 2 f 2 t -- 2 f 2 j -f ti f i t + fr i f r ) -Asin 2 9 k 59f rt \ . (4.7) 



We point out that the effective 6d metric that appears in this action is modified by the 
presence of worldvolume flux. Indeed, defining 

g^dxW = { 9lw - F w F1)dx"dx v = sin 2 9 k (ds 2 AdS2 + dil 2 2 ) (4.8) 



(4.7) becomes 



^ = /^^(-W^- \r u 9 pa f» P U + ^ 2 - 1 r{-5u 2 + 259 2 

sm 9 k J [2 2 sin 9 k 
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+ T 5 J d 6 Z(-4sm 2 e k 56f rt ). (4.9) 

The modification of the worldvolume metric by flux is familiar in noncommutative field 
theories obtained from B- fields, and is expected on general grounds [21]. By SO (3) sym- 
metry, 5u does not mix with the other modes; the mixed term between 56 and f rt agrees 
with the calculation in [20]- We also keep the other f ifi and fa components of the gauge 
field and we will explain shortly how to fix an adequate gauge. 

The next step is to perform the dimensional reduction, expanding the fluctuations in KK 
modes of 5 4 . The 6d scalars 5u, 56 and f rt are expanded in terms of scalar spherical 
harmonics, 

V%Y l (y) = -l(l + 3)Y l (y). (4.10) 

These modes transform in the /-th symmetric traceless representation of SO (5) (denoted 
by (0, 1) in terms of the Dynkin labels). On the other hand, the gauge field components Oj 
along 5 4 decouple from the 2d gauge field by choosing the gauge V*aj = 0. These modes 
are expanded in terms of transverse vector spherical harmonics, 

(V%-3)Y{(y) = -l(l + l)Yl(y), (4.11) 

which transform as the (2,1 — 2) of 5*0(5). Details of the dimensional reduction may be 



found in appendix A. 2 



4.3 Holographic dictionary for defect operators 

Let us summarize our results on the D5 fluctuations, and obtain the spectrum of defect 
operators. The fluctuations may be grouped in three sectors. 



• The (56, f r t) sector: these modes are mixed and, as we find in the Appendix and in |20j . 
their spectrum is 

m 2 J ^ ' ' = 1 ' 2 '- hl J 1 (4.12) 

1 l(Z + 3)(7 + 4) , I = 0,1,2,... 1/ + 4 V ; 



Their quantum numbers are (I + 4, 0; 0, 1) for I > 0, realizing the last state in (4.3), and 



(I, 0; 0, 1) for I > 1, corresponding to the first state in (4.3). In what follows, these two KK 
towers are denoted by (56, f r t)i and (56, frt)i respectively. 

• 5*0(3) triplet from D5 motion along du 2 + u 2 dQ 2 . The /-th KK mode has mass 

m 2 = 2 + 1(1 + 3) , I > 0, (4.13) 
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which is dual to an operator of dimension h = 1 + 2. These modes realize operators 
with quantum numbers (I + 2, 1; 0, 1) under SL(2, E) x SO (3) x SO (5) C OSp(A*\A). This 



corresponds to the third state in the supermultiplet (4.3), after identifying f = 1 — 1. These 
modes will be denoted by 5ui. 

• The KK modes from the internal gauge field a, give rise to SO (3) singlets with masses 



(see pUp ) 

mf = 1(1 + 1) , l>2. (4.14) 
The dual operators have quantum numbers (I + 1, 0; 2, 1 — 2), corresponding to the fourth 



state in (4.3). Below, we refer to these modes as (oj); 



Let us now establish the holographic dictionary for (chiral primary) defect operators. First, 
our analysis reveals that the lowest mode ((59, /rt)/=i i n the notation above) is dual to a 
marginal operator that transforms as an SO (5) vector. This is consistent with the ultra- 



short multiplet (4.2) of the superconformal algebra. Geometrically, this mode corresponds 
to the D5 scalar field that rotates the embedding SO (5) C SO (6). In the defect interaction 
L D x^'PaXi this scalar field gives rise to a fluctuation 5n a that is orthogonal to n a . Since 
an AdS field couples linearly to its dual operator, from this defect interaction we deduce 
that the operator dual to the massless scalar is X'P^Xi where n a <j)j_^ a = 0. 

At first sight, it may be a bit surprising that the only marginal operator that we find is an 
SO (5) vector and not a singlet. Indeed, naively we would expect the Kondo interaction 
x(Ao + n a (j) a )x, which is an SO (5) singlet, to give rise to a marginal deformation. However, 
the supersymmetry algebra fixes the coefficient of the Kondo coupling in terms of the 
ambient gauge coupling Qym- Thus, turning on x(A + na 4>a)X without changing g YM 
breaks supersymmetry and the dimension of the operator is not protected. Our results 
indicate that at large N this dimension diverges. 

The holographic dictionary may be found starting from the lowest weight mode (59, frt)^, 
which we just argued is dual to X^iXi an d then recognizing that the rest of the modes fill 



in a short multiplet of OSp(A*\A) -see (4.2) and (4.3). We summarize the dictionary for 



defect operators in ( 4.15| )J^| Only the bosonic fields are shown, and the operator content is 



schematic. The first two states realize the ultrashort multiplet, while the next four entries, 



defined only for I > 2, yield (4.3). 



5 Here a = 1,2,3 are the space coordinates of the four-dimensional spacetime, while the index i labels 
the S 4 coordinates (where the D5s are wrapped). Also, denotes the Af = 4 scalar fields that transform 
as an 50(5) vector, while the remaining scalar is n a (f> a . ((j)^ 1 ■ ■ -4>^) is a symmetric traceless product, 



and the symmetrization in (a^); corresponds to the vector harmonics of S , explained in [A. 2 
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D5 field 


defect operator 


SL(2,R) x 50(3) x 50(5) 


5ui =0 


o = x<f)±x 


(1,0; 0,1) 
(2,1; 0,0) 


(59, frt)^ 

(a>i)l 
(59, f rt )f ] 


o(^x(<f> i r---<f>i ) )x 

Q^^x(r t ny a <p [ l 1 ^ 2] --^ a l ) )x 
Q^ 1 ) ^x((n a D a <p a f^ ...ct>T 2) )x 


(l,0;0,l) 
(/ + 1,1;0,/-1) 
(/ + l,0;2,/-2) 
(/ + 2,0;0,/-2) 



(4.15) 



This concludes our analysis of chiral primary defect operators, that appear in the large TV 
spectrum. It is important to point out that, while the spectrum of fluctuations was derived 
in the probe approximation, the results are in fact exact. This is because the dimensions of 
chiral primaries are protected from quantum corrections. Next we will use these results to 
study the stability of the Kondo fixed point, and determine the thermodynamic properties 
induced by irrelevant deformations. 



5 Stability and defect thermodynamics 

So far we have studied the Kondo fixed point that arises from placing a fermionic impurity 
in M = 4 SYM , in a way that preserves half of the supersymmetries of the ambient theory. 
We determined the impurity entropy of the system and the spectrum of defect operators. 
The most interesting properties of the defect model are controlled by the lowest dimension 
operators allowed by symmetry that provide a correction to the scaling limit. In this 
section we use the results on the spectrum of operators to explore the stability of the fixed 
point and its basic thermodynamic observables. 



5.1 Stability of the Kondo fixed point 

The first point to address is whether the Kondo fixed point that we have found is stable; 
this would be the case if there are no relevant operators allowed by symmetries. It is 



useful to first review the situation in the real Kondo model (2.1), based on free electrons 
interacting with a spin. In the original one channel 377(2) spin model, the strongly coupled 
point J — > oo is stable, since by symmetry no relevant operators are allowed. In contrast, in 
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the multichannel case this infinite coupling point is unstable (the ground state degeneracy 
at strong coupling is larger than for J = 0). The model flows instead to an intermediate 
fixed point that is not described by Fermi-liquid theory. 



Let us ask the same question in our case. Given the results (4.15) on the spectrum of 
defect operators, the operator O = x4>±X °f lowest dimension is marginal but it transforms 
as a vector under 50(5). The lowest dimension singlet operator has then dimension 2, 
corresponding to O 2 . This is irrelevant in a + 1-dimensional theory. Thus, there are no 
relevant or marginal operators allowed by SO (3) x SO (5) and the Kondo fixed point is 
stable at large 't Hooft coupling. At small A, operators that are not chiral primaries may 
also have small dimension and should be taken into account. 

Notice that at large N and A, the stability of the fixed point is already guaranteed by 
the large 05p(4*|4) superconformal algebra. Indeed, in this limit the operators of lowest 
dimension arise from small representations, but this algebra does not have short multiplets 
with states of dimension h < 1 that are singlets under SO (3) x 5*0(5). This should be 
contrasted with known CFTs that contain light scalars and thus can have marginal or 
relevant deformations. The restrictions from the global symmetries are relaxed in systems 
with impurities separated along the spatial direction. This breaks 50(3), so operators 
with nonzero 50(3) spin are now allowed in studying the stability of the impurity fixed 



point. In this case the lowest allowed dimension is two, from the second state in (4.15). So 



we learn that the impurity fixed point in a system with broken 50(3) is still stable. 



5.2 Impurity specific heat and magnetic susceptibility 

Interesting observables in a defect model include quantities like the impurity free energy, 
specific heat and magnetic susceptibility. In the strict scaling limit, because the free energy 



of the defect will scale as infect ~ T (see (3.13)), the impurity specific heat 



a 



defect 



-T 



d 2 F r 



defect 



Q T 2 



(5.1) 



will vanish. 



The behavior of the defect susceptibility is more subtle and interesting. We define the 
susceptibility as coming from considering a background "magnetic field" B that couples to 
the unbroken 50(5) R-symmetry current via 



5 D / d A xA a J\ 



(5.2) 
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(with A the vector potential for B). The total susceptibility is then 

d 2 F 



Atotal 



(5.3) 

B=0 



ae 2 

and the impurity susceptibility is the piece of Xtotai that does not scale like a spatial 



volume, namely ^defect = Xtotai - Xambient- With our definition ( |5.2[ ) of magnetic field, the 
susceptibility is given by integrating the two-point function of SO (5) R-currents. In the 
original overscreened Kondo model the defect susceptibility vanishes at the fixed point, 
but this need not be the case when the ambient theory has interactions. For instance, 
in the models of [Hj, Xdefect = C/T with C a universal (irrational) number. For us 
Xdefect vanishes trivially in the probe limit; however, results with circular Wilson loops in 



6 6.4 suggest that it is nonzero after taking into account the defect backreaction, with a 



coefficient proportional to (3.14). Until we are ready to discuss the backreacted solution, 
we restrict the analysis of the susceptibility to the probe limit. In contrast, our results 
here for the specific heat are valid beyond the probe approximation. 

In the overscreened Kondo model, the defect specific heat and impurity entropy vanish in 
the strict scaling limit. It is because of this, as emphasized in e.g. [H[2], that the crucial role 
in governing the defect thermodynamics is played by the leading irrelevant operator that 
may be present in an RG flow to the critical point. In the overscreened model summarized 
this operator is a descendant Oq = J_i • (j) obtained by contracting the spin 



in 62.1 



current with the primary operator in the adjoint of the spin group. The dimension of Oq is 
h = 1 + N ^_ K , which becomes h = 3/2 for K = N. We stress again that for us Xdefect = 
only in the probe limit -below we find a nonzero susceptibility in the backreacted case, 
even in the absence of irrelevant deformations. 

We should therefore consider the results for Cdefect and Xdefect in the presence of such a 
perturbation. More precisely, one is generally interested in the leading operator consistent 
with the 5*0(3) x 5*0(5) global spin and flavor symmetries. Let us first discuss this for a 
general operator Oq that is a singlet under the global symmetries and has scaling dimension 
h® >2 (recall that there are no relevant or marginal singlet operators). In the dual field 
theory, we imagine deforming the defect action so that 

5defect ->■ 5 dc foct + J dt (\ O + h.c.) . (5.4) 

Here, Ao is a parameter of mass dimension 1 — ho < 0; we can use it to define a scale analo- 
gous to the Kondo temperature in [HI2], T K ~ Ag 1 ^ ■ The defect operators considered 
below have vanishing one point function so, as in the overscreened Kondo problem, the 
corrections to the specific heat and susceptibility are second order in Ao- Now, the basic 
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thermodynamic properties of the defect theory follow from dimensional analysis, recalling 
that the impurity specific heat is dimensionless and the susceptibility has dimensions of 
inverse length. At order Aq, the contributions then read 

r rp \ 2{ho-X) / T \ 2 ( fto_1 ) I 

Cdefect ~ 7fT~ ) Xdefect ~ JfT 7^ ■ (5-5) 



T K ' "™ \T K ) T 



Moving next to the specific form of the irrelevant perturbation, the analysis of KK fluctu- 
ations above shows that the lowest operator allowed by symmetries is the double-trace 

(xV-lX) ■ (x f 0iX) , (5-6) 

with a sum over the 5*0(5) indices of This operator has ho = 2, which curiously agrees 
with the normal Fermi liquid case (although in the Fermi liquid analysis the contribution 
was first order and not second order as here, so Cd e f ec t and Xdefect scale differently with 
temperature). On the other hand, the overscreened model relevant for us has a leading 
operator of fractional dimension ho = 3/2. This would give rise to a linear scaling Cdefect ~ 
T, as in the Fermi liquid case, although of course other properties of the system are sensitive 
to the non Fermi liquid nature of the fixed point. 

Large N correlators involving double trace operators are suppressed by additional powers 
of iV as compared to their single-trace counterparts. Thus, we may instead consider the 



leading single-trace singlet, which corresponds to the last state in (4.15). This operator 
is a descendant of = x(0±Vl 2 ' ) )x an d has dimension ho = 4. It is also interesting 
to envision the possibility that the two irrelevant operators are simultaneously present in 
the approach to the fixed point. This would produce a transition in the impurity critical 
exponents, at a parametrically small temperature controlled by a power of 1/N. 



6 Beyond the probe approximation: Kondo models and geomet- 
ric transitions 

Our analysis of the supersymmetric Kondo model has been restricted so far to the probe 
limit for the D5 branes. In the field theory side this corresponds to a quenched approxima- 
tion that ignores quantum effects from the defect fields, and in the gravity side it amounts 
to treating the D5 branes as probes in the D3 geometry. Moreoever, we found that even 
starting from this limit it is possible to derive results that are valid beyond the probe 
approximation -these include the spectrum of defect operators and the thermodynamic 
properties presented in £|5j 
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Nevertheless, it is important to go beyond the probe approximation and understand how 
the insertion of the impurity deforms the ambient gauge theory and its correlators. These 
effects become important for g\ M M > 1. What type of field theory do we expect after 
taking into account the backreaction from the impurity? In the boundary CFT approach, 
the impurity disappears completely and is replaced by appropriate boundary conditions on 
the ambient fields. In this section we study the backreacted solution for the supersymmetric 
model, and find a realization of this effect in terms of a geometric transition in the gravity 
description. The gravity solution can be trusted when all the cycles are large, which in 
particular requires g\ M M 3> 1. 



6.1 Basic properties of the backreacted solution 

Fortunately, the full solution including backreaction from the Wilson line defect is already 
known. The gravitational description of antisymmetric Wilson loops was first considered in 
[251 126] , and, impressively, the complete gravity solution was found by [7] . It is described 
by a "bubbling geometry" of the type introduced in [27J. Backreaction of the impurity 
can also be calculated in terms of a matrix model for the zero mode of A + n a <p a (when 
placing the theory on S 4 )^ Here we will summarize the gravity solution of [7J, reserving 
some comments on the matrix model approach to the last part of the section. 

Physically, taking into account backreaction replaces the D5 branes by 3-form field strengths 
F% and if 3, and the gravity solution does not include explicit source branes. This requires 
the appearance of new topologically nontrivial cycles ("bubbles") to support the fluxes. 
This is known as a geometric transition. Recalling the description of the AdS$ x S 5 limit 



in terms of (3.3) and since the full solution should still have SL(2,M) x 5*0(3) x 5*0(5) 



symmetry, the most general metric is of the form 

dn 2 

W = flds\ dS2 + HdVi\ + / 4 W + dY? , (6.1) 

where E is a Riemann surface with boundary (see below), and the radii fi are functions on 
E. Locally the space is of the form AdS2 x S 2 x S A x E, and the factors AdS2 X S 2 x S A are 
fibered over E. In the AdS§ x S 5 case, which we will revisit shortly, the Riemann surface 
is given by (u, 9), and c?E 2 = du 2 + d9 2 . 

In [7J it was found that the full solution is determined in terms of two real harmonic 



6 As recognized by [161 117] . matrix models provide powerful techniques to compute exactly certain 
quantities in M = 4 SYM. This connection was recently proved by [18] ; the matrix model for antisymmetric 
Wilson loops was developed in [22 j l8"l l28 l 129 ] 
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functions hi and /12 on E, given by the combinations 

hi = \e~*ftfl , hi = \e*flft , (6.2) 

where <fi is the dilaton. The /j and the dilaton are nonvanishing inside S, and the solution 
has a single asymptotic AdS§ x S 5 (corresponding to moving far away from the impurity). 
The boundary <9E is defined by the locus where either S 2 or S 4 shrink to zero. Since f 2 or 
/4 vanish on the boundary, h 2 = everywhere on <9£. Furthermore, in regular solutions f\ 
is nonvanishing on <9£; so hi = on segments on <9£ only when S 4 shrinks to zero. Given 
this structure, a topologically nontrivial 3-cycle is obtained by fibering S 2 over a one-cycle 
in £ that connects two boundary segments where f 2 = 0. Similarly, a homology 5-sphere 
can be constructed by fibering S 4 over a cycle between two boundary segments where 
/4 = 0. These are the cycles that support nonzero fluxes and represent the backreaction 



of the branes on the geometry. A concrete solution is studied in detail in £6.2, and the 
topology is illustrated in Fig. [2} 

Let us choose conformally flat coordinates on the Riemann surface, 

rfS 2 = 4p 2 dvdv . (6.3) 

It is also convenient to introduce the combinations 

W = dhidh 2 + dhidh 2 
V = dhidh 2 — dhidh 2 
Ni = 2hih 2 \dhi\ 2 - hjW 

N 2 = 2hih 2 \dh 2 \ 2 - h 2 2 W . (6.4) 
In terms of these quantities, the metric functions fi and p are given by [7] 



1/4 / \ 1/4 



Ni Ni J ■ \ V N 2 Z N : 

1/4 



'2 



U = 14,^ _™, . ,,5, 



iV! If I ' F ~ V h\h\ 

The dilaton (yf 2 = e 2 ^ is determined by 

e 2^ = _^ > . (6.6) 
Ni y ' 

In what follows we first formulate Ad S5 x S 5 in this language. Then explain how to choose 
hi and h 2 and consider and study the simplest case of one new 3-cycle. 
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6.1.1 AdS$ x S 5 warm-up 



In order to develop intuition about the boundary conditions that arise on <9£ and the form 
that hi and h% can take, it is instructive to first consider the AdS§ x S 5 case. In the general 
solution, this may be seen as the asymptotic description far away from the D5 insertions. 



The AdS$ x S 5 metric was written in the form (6.1) in £3.2 In this case the Riemann 
surface is a semi-cylinder < u < oo, — n/2 < 9 < n/2. Note here that at each point on 
the boundary of the S cylinder, the volume of either the S 2 or the S 14 vanishes. Along the 
locus 9 — ±| it is the S A that vanishes, while along the locus u = it is the S 2 . 

To make contact with the more general case, it is useful to define a new complex coordinate 
v via 

v = — i sinh(w + i9) . (6.7) 

The full boundary of the cylinder is now mapped to the line \mv = 0. For Rev < — 1 
and Rev > +1, f± vanishes; for —1 < Rev < +1, /2 = 0. The "central line" 9 = in 
the (it, 9) cylinder is mapped to the negative imaginary axis, and the whole (u, 9) cylinder 
is mapped to the lower half-plane in the v coordinate. The structure of the boundary is 
therefore a line along which loci with vanishing S 4 are separated by a locus with vanishing 
5* 2 . Drawing a curve into the interior connecting the distinct loci where the S A vanishes, 
one sweeps out a topologically non-trivial S* 5 . 

In order to verify that one has proper behavior at the asymptotics in more non-trivial cases 
with D5 branes, it is useful to map the asymptotic v — > oo region to a finite point. This 
can be accomplished by a final change of variables, defining: 

1 , . 

w - w = -— — . (6.8) 
1 + v 

Here, Wq is an arbitrary finite point, and w = Wq is where the AdS^ x S* 5 asymptotic region 
is mapped. So in the more general cases with D5 branes, even after backreaction, one will 
impose the condition that the solution asymptotes to AdS$ x S 5 as w — > wo. We will be 
more concrete about what this means, after developing a bit more technology. 

Note that in the new coordinates, if we choose wq to be real, we find that the region on 
the real axis where hi ^ in the f-variable, is mapped to the region (— oo, ei) on the real 
axis for the w- variable, with ei = Wq — 1/2. In fact, the Riemann surface £ can be viewed 
as the surface defined by the equation 

s(w) 2 = w — ei . (6.9) 
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We will henceforth always assume, in all cases, that wo is chosen to be real. 
In this notation, AdS§ x S 5 corresponds to 

hi = 2 coshw cos 9, h 2 = 2 smb. u cos 9 . (6.10) 
In terms of the v coordinates, this becomes 

hi = vl — v 2 + c.c., h 2 = i(v — v) . (6-11) 



Going through the final change of variables (6.8) obtains 



hi = \[2 — - yjw - ei + c.c, h 2 = i( h — — J . (6.12) 

W — Wo \ W — Wo W — Wo J 

Notice that the branch cuts are on the real axis, and that hi is precisely non-vanishing on 
the real axis when w > e\. 

In the remainder of the paper, we will switch back and forth between the v and w variables, 
depending on which is more convenient for a given discussion. 



6.2 General formalism for backreacted solutions 

This structure will generalise to cases including D5 stacks. The asymptotically AdS$ x 
S 5 nature of all of the relevant geometries means that the regions near v = ±oo will 
always have f± = 0; if there are p intervening regions where f 2 = 0, then there will be p 
topologically non-trivial S 5 s in the geometry, and p — 1 topologically non- trivial S 3 s. Not 
surprisingly, p — 1 then maps to the number of distinct D5 stacks (wrapping different polar 
angles along the S* 5 ) that we include, all sitting at the same point in IR 3 . We will focus 
on the case p = 2, so in the backreacted geometry, the Imt> = line will be split into 5 
regions: 3 segments with vanishing f±, and two interior regions with vanishing f 2 . The F% 
flux through the single non-trivial S* 3 will encode the number of D5s in the stack. 

Let us now describe in detail the case of a single stack of M D5 branes, namely when E is 



of genus 1. We see from (6.12) that in the simplest case (representing just an unperturbed 



AdS§ x S 5 geometry), one has 

dh 2 = i- -r . (6.13) 

{w - WqY 

That is, dh 2 has a double-pole at w = w - In fact, [7] found that this should always be the 
case for regular solutions that are asymptotically AdS§ x S 5 . So the only change from the 
backreaction of the D5 branes occurs in the harmonic function hi. 
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Now, as in §2, let us pick some angle 9 on the S 5 , and wrap M D5-branes there. Along the 
real w axis, we expect three segments where hi = 0, broken up by two finite regions where 
hi ^ 0. If we call the real coordinates where the transitions between hi — and hi ^ 
take place ei j2 ,3,4, then without loss of generality we can set e 4 = oo. The Riemann surface 
with boundary and the homology 3- and 5-cycles are shown in Fig. |2j 




s 5 



s ,r> 



s 3 



e 4 / e 3 e 2 \ t\ **., wg 

/ 4 = \ /a = _/7 4 = _ / 2 = / / 4 = 

Figure 2: Backreacted solution in the presence of a single D5 stack. The upper figure shows how 
the original cylinder corresponding to AdS^ x S 5 is deformed by the introduction of D5 branes, 
which create a new boundary segment with / 4 = 0. The lower figure shows the Riemann surface 
with boundary after the conformal transformations described in the text. The topologically 
nontrivial 3- and 5-cycles are shown in dotted blue lines. 

The defining equation of the genus one Riemann surface (where hi, h 2 are single- valued) 
is given by 

s 2 (w) = (w — ei)(w — 62) (w — e^) . (6.14) 

We can always redefine the w variable so that ei + e 2 + e% = 0, so there are two real 
parameters in the e^. These parameters will govern the integrated fluxes over the non- 
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trivial S 3 and S 5 s in the geometry; they are related to the number of D5 branes inserted, 
M, and the "shift from N" in the amount of F§ flux through the two S 5 s, which will be 
given by iV — k and k, with k determined in terms of the angle of the S 4 which the D5s 
wrap (see (3.6)). 



In the form (6.14), the surface is still non-compact. We can make it into a compact 
Riemann surface by going to Weierstrass form. In terms of the Weierstrass V function, we 
can define the z coordinate by the equation 

w = V(z). (6.15) 
Since the Weierstrass function identically satisfies 

3 

V\zf = 4 l[(V(z) - e t ) , (6.16) 

i=i 

we see that s(w) automatically satisfies s(w) = ^V'(z). 

Let us define the images in the z-plane of the special points in the w-plane via 

w = V(z ), a = V{uJi) . (6.17) 
Also, define the Weierstrass ( and a functions via 

C(z) = -V(z), C(z) = a ^. (6.18) 

Then the harmonic functions are given by [7j 

h = 2i U(z - zq) + ((z + z Q ) - 2^^^ - c.c^j (6.19) 

and 

h 2 = (C(* + zq) - ((z - Zq) - c.c.) . (6.20) 

This is the desired final expression for h\ and which provides an explicit solution for 
the single D5-stack case. Generalizations to more stacks can be similarly worked out in 
terms of hyperelliptic Riemann surfaces and harmonic functions on them. 

The number of D5-branes and the 9 angle of the S 4 in S 5 which they wrap, should be 
indicated by the amounts of three-form and five-form flux threading the various topologi- 
cally non-trivial spheres in the geometry. In terms of the holomorphic and antiholomorphic 
parts 

h 1 = A + A, h 2 = B + B, (6.21) 
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the expressions governing the fluxes are as follows [28] : 



J F 3 = 47rVM = 87r / (idA + c.c.) (6.22) 



(-2 

and 

"e-21-l 



r re-21-i 

/ F 5 = 4n 4 (a') 2 N I = 8n 2 (AdB - BOA + c.c 

Js 5 r Je 2I 



(6.23) 



In (6.23), the two different 5-spheres are labeled by / = 1,2, and the integration contour 
should actually be taken just a bit below the real line segment running from e^i to e2/-i - 
i.e., with small negative imaginary component. 



6.3 Matrix model approach 

We just described the solution very explicitly in terms of elliptic functions and geometric 
transitions where new topologically nontrivial cycles appear. At this stage, it is natural to 
ask how all these properties are understood in field theory language. The simplest answer 
is given in terms of a matrix model that keeps the zero mode of A + <ft when the theory is 
placed on S 4 . Intuitively, the distribution of eigenvalues in the matrix model will describe 
the homology cycles, and the Riemann surface will correspond to the spectral parameter. 
Let us briefly explain this approach, following [H [28j 129] . 

The claim is as follows. Let M in this section denote an NxN Hermitian matrix associated 
to the zero mode of A + n a (j) a ^ The expectation value of a gauge-invariant function f(M) 
-i.e. fiVMV^ 1 ) = f(M) for a unitary matrix V- is given by 

(/(M)) = iy dM f(M) exp ^-^tr[M 2 ]) , (6.24) 

where the partition function Z is defined by 

Z = j dM exp ^-^tr[M 2 ]^ . (6.25) 

Here the measure dM denotes Yl dMij subject to the constraint of hermiticity. Then, the 
claim is that the BPS Wilson loop (whose action is identical to that discussed in §2.2), 
in the R representation of the gauge group, can be calculated in the matrix model by 
computing (tr R [e M }). 



7 This is not to be confused with a number of D5-branes, which was denoted by M in earlier sections. 
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The partition function reduces to an integral of the exponential multipled by a familiar 
Vandermonde determinant encoding the eigenvalue repulsion of the eigenvalues rrii of M: 



N / N 



Z ~ J \\drrij Tk<i<j<N (m, - mrfexp ^--^5Z m ij " (Q.2Q) 

As usual, the eigenvalue repulsion can be exponentiated into the "action" in the matrix 
integral, yielding a logarithmic potential: 

N 

Z~ Yl dm i exp(-S(m)) (6.27) 



i=l 

TV 

A" 



S'(m) = — ^- — 2 log|mj — m,j | . (6.28) 

j=l i<j 

Consider now the insertion of the Wilson loop in the k-th antisymmetric representation; 
after diagonalizing M, the trace becomes 

trA fc [e M ] = ex ^ m n + - + m h) ■ (6-29) 

I<jl<j2.<..<jfc</V 

As a result, we can also absorb the operator insertion into a new effective action and do a 
purely exponential integral: 

(trAje A/ ]> ~ ( {N W k)]k] ) J A dm i eMSk(M)) (6.30) 



where 

S k (M) = S(M) -J2 m >- ( 6 - 31 ' 



=i 



Introducing the resolvent 

w (*)= tr 7=M' ( 6 - 32 ) 

one can show that the equation of motion for m; relates u to the harmonic functions of 
the bubbling solution via [281 I2H] 

fti oc 2(z — z) — (u) — lo) , h 2 oc z — z (6.33) 

once the spectral parameter z is identified with the Riemann surface coordinate v. In this 
way, the geometric transition can be seen directly in the field theory]^] 

8 In fact, it is easy to evaluate the integral by saddle point in the probe approximation, and solve for 
Sk [H]. The result agrees with the probe calculation in the gravity solution. 
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6.4 Comments on correlation functions 



Having described how to calculate quantum corrections from backreaction of the defect, 
the next step is to compute bulk correlation functions. The simplest correlators are one- 
point functions which, in the gravity side, follow from expanding the background solution 
around the asymptotic AdS^ x S 5 region, and extracting the leading corrections. Conformal 
invariance restricts the expectation value of a chiral primary Oa of dimension A to be of 
the form ^ 

(0*(x)) w = (6.34) 

where the subindex W denotes the correlator in the presence of the Wilson line insertion, 
I is the distance between the line and the operator position, and Sa is a constant that 
depends on A and the representation of the Wilson loop. All the nontrivial information 
is thus in Sa, which can be calculated either using the gravity solution or in terms of the 
matrix model. One-point functions of chiral primaries in the presence of a Wilson loop in 
the antisymmetric representation were calculated in [29J. 

Higher correlators are harder to compute and, to our knowledge, have not been obtained. It 
would be quite useful to fill in this gap, as these correlators contain important information 
about the response of M = 4 SYM to a nonlocal observable such as the Wilson loop. Let 
us outline the first steps in this program, which are also necessary for calculating trans- 
port properties away from the probe limit. We should start by understanding correlation 
functions of the gauge invariant chiral primary of lowest dimension, which is 

2 = ^ tr (W - ^ a V> a ) (6.35) 

where here a, b are 50(6) indices, and the normalization is such that in the planar limit 
(C 2 (a;)(9 2 (y)) = l/\x — y\ 4 . This dimension 2 operator is very important, because its 
descendants include the R-symmetry current (J^) a b and energy momentum tensor T^ u . 
Therefore, the correlation functions for (9 2 determine, via super conformal transformations, 
the basic properties of the correlators of and T^ u [30] . The one point function coefficient 



in (6.34) is given by [29J 

\l/2 

S 2 = — sin 3 4, (6.36) 

67T 



which is proportional to the g- function (see (3.14)) 



Approximate correlation functions can be obtained if we consider a circular Wilson loop 
(instead of a line operator), and expand it in terms of local operators [31] 



^ circ ie(a) = (W circle ) £a (27ra) A O A (0) (6.37) 
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with a the radius of the circle. This is valid for insertions at points far away from the circle. 
In this way, correlation functions with a Wilson loop insertion are approximated by higher 



point functions in the unperturbed TV = 4 theory. Since a comparison between (6.37) and 



(6.34) reveals that for chiral primaries £a oc Ha, in this approximation the correlators are 



determined by the one point functions that were discussed before. 

As an application, consider the correlator {J^(x)J u (x'))w for SO (6) R-currents in the 



presence of a circular Wilson loop. Expanding the Wilson loop as in (6.37), the lowest 
dimension contribution is proportional to a 2 H2 ^J^JyO^. The M = 4 correlator that we 
are left with can be obtained by applying superconformal transformations to the 3-point 
function of Oi- This suggests that also the current-current correlator in the presence of a 



Polyakov loop will be non-zero, with Xdefect ~ H 2 /T which, given (6.36), is proportional to 
the irrational ^-function. This is reminiscent of the behavior found in [14J . On this note we 
end our anaysis, and hope to come back to the determination of correlators and transport 
properties in a future work. 



7 Discussion 



We have seen in this paper that the most basic physical questions about the maximally 
supersymmetric Kondo problem can be answered by standard holographic techniques. The 
operator content of the theory is determined by harmonic analysis on the geometry wrapped 
by the probe branes which represent the impurity in the gravity dual, and the backreaction 
of the impurity on the bulk liquid is encoded in a geometric transition. Several of the prop- 
erties of our model show striking resemblance to results on multi-channel Kondo problems 
in the field theory literature, but there are important differences, stemming from the power 
of the superconformal algebra which controls the IR fixed point. For instance, although 
some defect thermodynamic properties are most analogous to those of overscreened multi- 
channel Kondo models, the leading irrelevant operator at the fixed point (which governs 
some of the most interesting physics of the defect) is dimension two, as in the original 
single-channel Kondo problem. 

It is certainly of interest to extend the results of this paper to systems with less symmetry. 
Clear directions for future research include: 

• While several features of the model seem analogous to those of multi-channel Kondo mod- 
els studied in the condensed matter literature, some important differences arise. Notably, 
instead of obtaining non-Fermi liquid scaling of the leading irrelevant operator, we obtain 
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a dimension 2 operator, precisely the same as in the normal Fermi-liquid case (neverthe- 
less, the contribution is second order in the coupling). It seems likely that this is because 
the superconformal algebra heavily constrains the spectrum of operators, resulting in the 
leading irrelevant A = 2 operator we described in §5; certainly, the itinerant sector here 
has little resemblance to a single-channel bulk Fermi liquid! In similar models of impurities 
interacting with less symmetric CFTs than M = 4 supersymmetric Yang- Mills, we antici- 
pate that the results will be less constrained and that the coincidence in the behavior of 
the leading irrelevant operator will disappear. 

• We have studied here the interactions between the defect and the ambient theory at 
vanishing chemical potential. A more realistic analogue of the Kondo problem would be 
to consider the M = 4 theory at finite /i for the 5*0(6) R-charge; even more similarity to 
the real problem could be obtained by using the ideas in [32] to model a bulk Fermi liquid. 
Needless to say, finding the finite temperature solutions could also be of interest. 

• It would be very useful to study correlators of Af = 4 SYM in the presence of a Wilson 
loop. To our knowledge, this has not been done yet. Furthermore, with the backreacted 
gravity solution in hand, we can now ask how the defect modifies bulk transport. For 
instance, we can study the bulk resistivity by computing current-current correlators in the 
modified background holographically. Such studies, either in this microscopic model or in 
"bottom-up" constructions of similar models, might allow one to find interesting general- 
izations of the original "resistance minimum" that motivated the tremendous theoretical 
interest in the Kondo model to begin with. 

• Our results motivate further study of overscreened models with equal number N = K 
of spins and channels. In this case the ambient electrons become N x N matrices, and we 
described powerful matrix model techniques to solve this system. Studying more general 
large N matrix models for the Kondo problem could add a valuable tool to the techniques 
that already exist in the literature. Even at the level of the M = 4 SYM, the matrix model 
description should be generalized to include defect operator insertions. 

• Finally, and most ambitiously, one would like to push for a much more complete un- 
derstanding of lattices of such defects, including their backreaction on the bulk and their 
interactions with one another. The results in (H [10] demonstrate that the simplest su- 
persymmetric lattice models can give rise to interesting non-Fermi liquid physics in the 
semi-holographic approach [33J. But the analysis in (TU] also suggests that backreaction 
at some power-law suppressed scale in 1/N will destabilize the non- Fermi liquid (as it de- 
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stroys the AdS2 geometries which are at the root of the non- Fermi liquid behaviour). Our 
analysis here sheds interesting light on those results. A priori, one might have thought 
that inter-defect interactions would lead to multiple defects "destabilising" each other's 
Kondo fixed point. In RG language, the lattice model breaks the £0(3) symmetry (but 
not necessarily the 5*0(5) symmetry) enjoyed by the single-defect fixed point. Hence, 
defect-defect interactions could a priori induce relevant S'0(3)-breaking operators in the 
fixed point Lagrangian, destabilizing the fixed point. However, a surprising result of our 
analysis in §4 is that even with SO (3) breaking, there are no relevant defect operators at 
the fixed point. This is not inconsistent with the discussion in [TO]; in fact, it is not inter- 
defect interactions, but defect interactions with the bulk that must destabilize the bulk 
phase in the supersymmetric Kondo lattice. This is plausible because the N = 4 theory 
has a marginal operator that respects all relevant symmetries; in the dual gravity, it maps 
to the dilaton, and gravitational backreaction can source the dilaton, driving one out of 
the gravity regime entirely. This does suggest, however, that if one can construct similar 
defect models where the ambient theory of the itinerant modes does not have marginal 
operators respecting all of the relevant symmetries, one may be able to find lattice models 
which remain stable, with AdS% regions on the lattice of impurity branes, order by order 
in 1 /N. 
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A Spectrum of fluctuations of the D5s 



In this Appendix we calculate the spectrum of fluctuations for the 6d theory on the D5 
branes, placed in the background of the color D3 branes. We will first compute the 



quadratic 6d action from (3.7), and in a second step perform the dimensional reduction 
on 5 4 . The relevant quantities are the background 6d metric (3.8), the worldvolume flux 



(3.13), the 4-form flux (3.4), and the embedding (3.5). 



The bosonic fields on the 6d worldvolume are a gauge field and the scalar fluctuations 
around the embedding (3.5 ). These are 9 = 9 k +S9, u = 5u, and two angle fluctuations along 
the 5 2 . Near u = we have du 2 + sinh 2 u dVt 2 ~ ds 2 ^, so Su and the two angle fluctuations 
combine into a 3- vector (p m . These transform as a fundamental of 5*0(3), which becomes 
the R-symmetry of the 6d theory. In the original 4d spacetime, ip m is simply the fluctuation 
in the 3d spatial position of the defect. For the purpose of computing the spectrum, it is 
enough to study 5u, since the masses of the triplet are then the same based on the 5*0(3) 
symmetry. 



A.l Six-dimensional quadratic action 

Let's start by expanding the determinant term of the DBI action in the fluctuations. Define 
M^ u = G^y + F^ u . Then the background value is 



/ - 



1 cos t 



M 



Olx4 \ 
Olx4 

V 4x i 4x i sin 2 9 k g^j 



COS tffc 1 



(A.l) 



where the zero entries are blocks of the indicated size. The expansion of the first term of 
the DBI action gives 

1 + -ti(M~ l 5M) + -(tr(M -1 5M)) 2 - -tr (M^dM) 2 
2 8 4 

(A.2) 



- det (M + SM) 



det M 



The metric fluctuation reads 

6(ds 2 ) = 6u 2 ds 2 AdS2 + (sin 29 k 59 + cos 29 k 59 2 )dtt\ + (ddu) 2 + (d59) 2 , (A.3) 

and for the field strength we write 

„ „ dt Adr , 

F = cos9 k =— + /, (A.4) 
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where / is a 2-form along 6d space. Thus, 

SM^ = (sin 26 k 56 + cos 26 k 56 2 )g^ + 5u 2 + d^(5u)d u (5u) + 8^56)3,(56) + 



Raising an index we obtain 
1 
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sin 2 6 k 



\ 04x1 04x1 ,. ,, 

I r 2 [-U u + cos 6 k U rt ] r 2 [cos 6 k U, 

„2 r n tt , tt 1 „2 rrr 



o o o lx4 

\0axi 4x i 5) ) 
,. T - U rt ) r 2 [cos 6 k U ri - U ti ]\ 



r 2 [- cos 6 k U tt + U rt ] r 2 [U rr - cos 6 k U rt ] r 2 [U ri - cos 6 k U u 

w 



\ u\ U{ u\ 

Here U^ v = d ^(511)8,(511) + d^(56)d u (56), and upper j indices are raised with g( s *\ 
Plugging (A. 6) into (A.2) obtains, up to quadratic order, 



V- det M = 



sm°6 k 2 
-5u — 



2 r 



:,,3 n 

~ [r 2 {-(d t <P a ) 2 + (d r <P a ) 2 ) + d^aO^a] + 



- I s in 3 4cos^(r 2 / rt + 4sin^^) + ^ [ r ^f^- ff. + r 2 (-f ti f i t + f ri f r )] + 

- -^(1 + 2 cos26 k ) sin 3 6 k 56 2 -4sin 2 ^ cos 2 6 k 56 f rt . (A.6) 
Here the kinetic term for a stands for the sum over the kinetic terms for 5u and 56. 



Next we consider the CS term; up to quadratic order, it gives 
FAC 4 = <f£yfg^ 



4^^ sin 4 6 k 56 - (h k - sin 20* + \ sinA6 k ) f rt 



4 sin 4 6 k 56 f rt + — cos 2 6 k sin 3 6 k 56 2 



(A.7) 



Finally, combining (A.6) and (A.7) obtains 



sin 3 6 k 1 



S {2) = T 5 / d^^) { [r 2 (d t <p a ) 2 - (d r <p a ) 2 - tfMfa - 25u 2 + A56 2 ] 



+ 



sin 6 



k / J f2 



) I r frt " 2^ # - ftift + frifr ' ~ 4 8 ^ ^ fH 



(A.8) 



This is our final 6d expression. (The term linear in 56 cancels, while the linear term in f rt 
is a total derivative.) 
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A. 2 Dimensional reduction 



SO (3) triplet 

We first consider the SO (3) triplet; we just found that it decouples from the rest of the 
modes, and this is consistent with the 5*0(3) symmetry. It is enough to consider the 
magnitude of the 3-vector, Su; it is expanded in spherical harmonics 



8u = ^5ui{a)Y\y) , V%Y l = -/(/ + 3)Y l . 



(A.9) 



Here a are coordinates along AdS2, y are coordinates on S" 4 , and I is a multi- index (see 



below). The action follows from (A. 8): 

S D T 5 sin 3 # \ J d2(J ^ (r 2 {d t 5u lm f - r 2 (d r Suim) 2 ~ W + 3) + 2) 8uf m ) . (A.10) 

lm 

Then the mass of the Z-th KK mode of the SO (3) triplet and the dual dimension are 

mf = 2 + 1(1 + 3) , hi = 1 + 2. (A.ll) 



Internal gauge field 

The 6d gauge field gives a M = (a t ,a r ,ai), where a« are along S 4 . We now consider these 
internal components, that reduce to scalars along AdS2- In order to decouple them from 
(a t ,a r ), we choose transverse gauge 

V*ai = 0. (A.12) 

Then, for instance, 

fuft -> -d t aid t a l - dia t d' l a t . 

The cross-term dia t d t a % has disappeared after integration by parts and use of the gauge 
fixing condition. 

Then the action for a, simply becomes 

S D r 5 sin0 fc J d Q i l - {[d tai ) 2 - (d r a z ) 2 + ^a, (gfaVfa) - ^ 4) ) aA (A.13) 

where i,j indices are raised and lowered with gis A )t an d the curvature of the internal sphere 
comes from integrating by parts the term fijf 1 -' . For a unit-radius S n , R 1 ^ = (n — l)g^ . 
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The KK reduction is performed by expanding etj in vector spherical harmonics on S 4 . 
Let us review the construction of vector harmonics for S n ^\ Start from the fundamental 
representation Y A of SO(n + 1), A = 1, . . . ,n + 1, which satisfies 

V 2 Y A = —nY A . 

Spherical harmonics are gien by symmetric traceless products: 

Y A 1 ...A l = Y ( Al _ _ _ Y A t ) ? V 2yA 1 ...A ! = + n _ yyAL.-A, _ 

These are the harmonics that were denoted by Y l above. 
Next we introduce 

Y A = d l Y A . (A. 14) 

Some useful properties are 

Y,Y i A Y? = -Y A Y B + 8 AB , J2( YA ) 2 = 1 > ViY 3 A = -5ijY A . (A.15) 

i A 

Notice that Y A is an eigenvector of the laplacian, but it is not transverse. 

The transverse vector harmonics are constructed as 

Y M... Al = y [a 1y{ a 2 ) _ _ _ yAl ) ^ i> 2 , (A.16) 

where the first two indices are antisymmetrized and the remaining / — 1 ones are sym- 
metrized. This gives a representation of SO(n + 1) with / — 1 boxes in the first row and 
one in the second -this will be denoted by (2,/ — 2). The calculation of the eigenvalues 
gives 

(V|n - (n - 1)) Y Al - Al = -1(1 + n- 3)Y Al ~ A < (A.17) 

Now we are ready to perform the dimensional reduction, specializing the above results to 
n = 4. Expanding 

a i (a,y) = J2*i(v)Yl(y) (A.18) 
(where / is a multi- index standing for A±, . . . , A{) we find 

SD±T 5 sm6 k J2 I d^^{r 2 (d t a{) 2 -r\d r a{) 2 - 1(1 + 1)^) , ^>2. (A.19) 

So the (2, / — 2) vector harmonic has a mass m 2 = 1(1 + 1) and is dual to an operator of 
dimension hi = I + 1. 

9 We thank N. Bobev for explanations on vector harmonics. 
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(59, frt) sector 



In this sector the action is 



S D T 5 J d^Vgl jsin 3 9 k {^(d t 89f- l -(d r 69f - ^(di69) 2 + ^69 2 ^j 
+ S * n ® k r 2 ( (d r a t ) 2 + (d t a r ) 2 — 2d r a t d t a r + ^-(d i a t d' t a t — dia r d l a r 



- 4 sin 2 9 k 59(d r a t - d t a r ) } . (A.20) 

We still have gauge redundancies along the 2d space, and we use them to set a Q = 0, after 
imposing Gauss' law: 

d r (r 2 d t a r + 4 sin 9 k 69) = 0. (A.21) 
Then the equations of motion simplify to 

2 

r 2 (-d 2 t 69 + d 2 .69) + V%69 + 459 + A^—d t a r = 

sin 9 k 

-d t (r 2 d t a r + Asm9 k 59) + V 2 s ±a r = 0. (A.22) 
Defining the combination 

n = -^--r 2 d t a r + A59 (A.23) 
smt^c 

and expanding in spherical harmonics, obtains 

d r T]l = 

r 2 (-d 2 59i + d 2 r 89i) - 1(1 + 3)56i - 1266i +4^ = 

-r 2 d 2 m - 1(1 + 3)( Vl - U9{) = • (A.24) 

The last equation here follows by taking one time derivative on the a r equation above. 

For the I = mode, this system gives d t f]o = d r r]o = 0, so that r] = const. Then the 
equation for 59 gives m 2 , = 12, which is dual to an operator of dimension ho = 4. For the 



higher I > 1 modes, diagonalizing the mass matrix in (A.24) obtains 



in agreement with [20]. This ends our analysis of the D5 fluctuations. 
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